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I. INTRODUCTION 



Because of asymptotic freedom, one thinks that the structure and properties of high- 
temperature quantum chromodynamics (QCD) should be accessible by perturbative means, 
however these are organized [l| . Originally, the standard loop-expansion ran into diffi- 
culties when attempting to describe the properties of slow-moving particles: gauge invariance 
was lost , a consequence of the fact that the standard loop-expansion did not reflect any- 
more a systematic expansion in powers of the small coupling constant. A re-organization 
of lowest-order terms in the correlation functions was therefore necessary and indeed per- 
formed: the so-called hard thermal loops (HTL) were systematically extracted from one-loop 
order diagrams and added to the lowest-order quantities HHQ- It should 

also be mentioned that the need for resummation of standard loop-expansion is not peculiar 
to non-abelian theories but is also necessary in the Abelian case, for quantum electrody- 
namics (QED) [12] and scalar QED as well as for the simpler <£> 4 -theory Q|. Thus are 
formed dressed propagators and vertices and a convenient description of the slow-moving 
quasi-particles becomes possible. Next-to-leading order quantities, in this new framework, 
are obtained via one-loop diagrams involving these dressed propagators and vertices. In this 
new context, some of the old difficulties are cured but others persist, most notably sensitivity 
to the infrared and the light-cone, essentially due to the fact that static (chromo) magnetic 
fields remain unshielded even after the HTL reorganization. In QCD, these are thought 
of becoming screened at the so-called magnetic scale, which is believed to manifest itself 
non-perturbatively in the dispersion relations at next-to-leading order. Because of this, per- 
turbation is believed to break down at some order, which depends on the quantity under 

EH 



consideration 



121 ■ 



Still, useful information can be extracted from the perturbative regime. Indeed, using 
this method, an early explicit calculation of the zero-momentum transverse gluon damping 
rate shows that it is finite and positive jisl l . A similar computation for the quarks has also 
been carried out independently in Q andQ and it yields a finite positive number too. As 
for the calculation of static quantities such as free energies and screening masses, these are 
most conveniently carried out in the imaginary-time formalism since analytic continuation 
to real time is not necessary E] anymore. For such quantities, there exists a simplified 
resummation scheme Eh based on the fact that in Euclidean formalism, the momentum 



2 



can be soft only for the zero mode in the Matsubara sum. Hence dressing is necessary only 
for the static propagators, and calculations based on this are relatively simpler to carry out 
than those of the full resummation, necessary when considering dynamic quantities. The 
'reduced' approach has been used in many instances like in the study of phase transitions in 
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2J as well 



gauge theories |22| , the calculation of the free energy in yj 4 -theory and QCD 
as in the computation of the electric screening mass in QED and scalar QED 21]. Other 
important quantities in this regard are the next-to-leading order Debye screening length 2^ 
and next-to-leading order correction to the gluonic plasma frequency 

But as mentioned, difficulties still persist after HTL reorganization. For example, log- 
arithmic sensitivity is encountered in an early estimation of the damping rate for a heavy 



fermion |8j], and more generally for fast-moving particles 27], |28|, |29|, |30| ■ Also, an estima- 



ion of the damping of moving quasi-particles is logarithmically dependent on the coupling 
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. The infrared problem was also emphasized when attempting to estimate the damp- 
ing rates of non-moving longitudinal gluons 33, 34] and slow- moving quarks 35]. 

What is of concern here is, partly, the analytic properties of QCD at finite temperature 
in the infrared. In order to investigate some of these properties, the context of scalar QED is 
useful because it offers a simpler setting in which calculations are facilitated by the fact that 
there are no hard thermal loops in the vertices and those related to the scalar self-energy are 
momentum independent kJ. This situation allows sometimes carrying out almost complete 
momentum-dependent calculations. Of particular interest to us is an issue raised by the 
works in order to carry forward analytically, one had to expand in powers of 

the external momentum p, considered ultrasoft to ensure the validity of the expansion, to 
allow one to perform intermediary solid-angle integrals which would have been, otherwise, 
intractable analytically. The particular feature of this expansion is that it was done in 
the imaginary-time formalism, before the Matsubara sum was performed and the analytic 
continuation to real energies taken. The subsequent appearance of infrared sensitivity in 
the final coefficients may then be linked to this 'early expansion'. 

To tackle this specific issue in more depth, we embark in this work on the investigation 
of the next-to-leading order dispersion relation for scalars with ultrasoft momentum p in the 
context of next-to-leading order HTL perturbation of scalar QED. We will examine both the 
damping rate 7 S (p) and the energy oj s (p) . Note that the damping rate 7 S (p) has already 
been investigated in and found, to second order in p, logarithmically infrared sensitive. 
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In this work, the expansion of 7 S (p) is pushed to fourth order in p and that of u> s (p) to second 
order. What is interesting is that we perform the expansions by two different methods: (i) 
We delay the momentum expression until after the Matsubara sum is performed and the 
analytic continuation to real energies taken. As for the damping rate 7 S (p), we will find the 
coefficient of second order logarithmically sensitive to the infrared cut-off 77 as in 3y] whereas 
the coefficient of fourth order will behave like I/77 2 (odd-order coefficients cancel). As for 
the energy uj s (p), no infrared sensitivity will appear, (ii) We perform the early- momentum 
expansion and carry on with steps similar to pa . [mJ . We find exactly the same results as 
those of the previous method. This must alleviate some pressure on the early-momentum 
expansion method: it may not be responsible for the appearance of infrared sensitivity. In 
the context of QCD, there is simply little hope of obtaining compact analytic expressions, 
and expansions such as the early-momentum one are sometimes necessary. 

This article is organized as follows. After this introduction, the second section recapit- 
ulates the hard thermal loops and the dressing of the propagators. Section three discusses 
the calculation of the scalar damping rate j s (p) by the two methods and section four the 
calculation of the scalar energy u s (p). Section five is devoted to the discussion of the results 
and includes some concluding remarks. 



II. DRESSING THE PROPAGATORS AND THE DISPERSION RELATION 

Let us first recapitulate the results giving the propagators and vertices to be used in the 
forthcoming next-to-leading order calculations. We use the imaginary-time formalism in 
which the four- momentum is P M = (po, p) such that P 2 = Po + p 2 with the scalar Matsubara 
frequency po = 2imT ', n an integer and T the temperature. After the evaluation of all 
frequency sums, po is analytically continued to the real external energy u using the analytic 
continuation p = —iu + + . 



A. Dressed propagators 

To leading order, the self-energies are obtained from undressed one-loop diagrams. For 
soft external momenta, that is for to and p of order eT where e is the coupling constant, the 
dominant contributions to these diagrams come from loop momenta of the order of T, the 
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hard scale. These contributions are called hard thermal loops. For the scalar self-energy, 
the hard thermal loop is given by 

5E(P) = -m 2 s , (2.1) 

with m s = eT/2, the scalar thermal mass. Note that is momentum independent. When 
dressed with the hard thermal loop, the leading-order propagator for the scalar becomes: 

(P) = • (2-2) 

s 



P 2 + 



m 



The hard thermal loop SH^ in the photon self-energy is given in |3i| and can be expressed 
in terms of two independent scalar functions 511; (K) and 8Tlt(K) given by: 



ike 

k J ' 5' 



5U l (K)=3m 2 p Q 1 [^); 5U t (K) = -m\ 



(2.3) 



where m p = eT/3 is the photon thermal mass and the Q^s are Legendre functions of the 
second kind. In the strict Coulomb gauge, the components of the dressed photon propagator 
*A fll/ (K) are as follows: 

*A 00 (K) = *A t (K) , *A 0l (K) = 0; 

*A;, (K) = (8^ - kk) * A t ( K ) > (2-4) 

where *A; and *A t are the propagators for the longitudinal and transverse photons respec- 
tively. They have the following expressions: 

' A >w= v-L l{K y " A ' {K)= w^km- (2 - 5) 

Now one peculiarity of scalar QED is that the vertices remain undressed, i.e., unaffected 

n 

by the hard thermal loops 13]. The vertex with one photon and two scalar external lines 
(Q incoming, P outgoing) is: 

T»(P,Q) = -e(P + Qr , (2.6) 
and the vertex between two photons and two scalars is: 

r""(P,Q) = 2e 2 5> MV . (2.7) 
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B. Dispersion relation 



The scalar damping rate 7 S (p) and energy uj s (p) are obtained from the scalar complex 
energy uj which satisfies the following full scalar dispersion relation js^J: 

cu 2 =p 2 -£(^,p), (2.8) 

where £(u;,p) is the full scalar self-energy Taking account of the next-to-leading order 
terms in e and expanding everything around the leading order result u s0 (p) (squared) for 
fixed p, we obtain: 



uj 2 



p 2 - (uj s0 ,p) - *£ (u s0 ,p) - (uj 2 - uj 2 s0 ) <9 X . 2 <5£ (x,p)\ x=ulaQ + O (e 4 T 2 ) . (2.9) 



In this relation, uj s0 (p) = ^Jm 2 s + p 2 , pole of *A S (P) from ()2.2j) and solution to (|2.8j) to 
lowest order. is nothing but the hard thermal loop given in (j2.1j) and *£ is the next-to- 
leading order contribution. Using (j2.1J) . we can rewrite (|2.9|) as: 



uj = UJ 



sO 



- *Z(uj s0 ,p) + O(e 4 T 2 ). (2.10) 



The full energy uj (p) is in general complex. If we denote its real part (the scalar energy) 
by uj s (p), then we have: 

Re5*Z{uj s0 ,p) , 3 . 
u; s (p)=w s0 = YO\eT). (2.11) 

ZUJ S Q 

The prefix 5 before *£ in f!2.11|) indicates that the known leading hard-thermal- loop contri- 
bution has to be subtracted for ultraviolet convergence jl^. The damping rate for scalars 
is defined by y s (p) = —Imuj (p). It is e-times smaller than uj s0 (p), and so we have to lowest 
order j^: 

7 S (p) = ^- Im *£ (u; s0 , p) + O (e 3 T) . (2.12) 

We see that determining y s (p) to lowest order in e and oj s (p) to next-to-leading order 
amounts to calculating the imaginary and real parts of the next-to-leading order scalar self- 
energy *£. From now on, we assume the scalar momentum p ultrasoft, i.e., of the order of 
e 2 T. Also, we will take m s = 1 in the sequel, to ease the notation. We will reintroduce it 
back in the final expressions. 



6 



III. SCALAR DAMPING RATE 



We start by ca 
two contributions 



culating the damping rate 7 S (p). The self-energy *E (P) is the sum of 



*E(P) = *E 1 {P) + *E 2 (P), (3.1) 



where the contribution *Y>i is a one- loop dressed diagram with two one-photon-two-scalar 
vertices given in (|2.6|) : 

*Ei (P) = Tr soft [T" (P, Q) *A S (Q) I* (Q, P) (if)] , (3.2) 

and the contribution *E 2 is the dressed tadpole diagram with a two-photon-two-scalar vertex 
given in ()2.7|) : 

*E 2 (P) = Tr soft [r^ (if, P) * A^ (if)] . (3.3) 

In the above two relations, if is the internal photon loop-momentum, Q = P — if and 
Tr = T^2 ka j d 3 k /(27r) 3 with fc = 2nirT. The subscript 'soft' indicates that only soft 
momenta are allowed in the integral. To take account of potential infrared sensitivity, we 
will introduce an infrared cutoff rj in the ^-integration. 

It turns out that *E 2 (P) is real. This is because there is only one dressed propagator 
involved in its expression ()3.3j) and the vertex there is undressed. It will also turn out to be 
p-independent, see (|4.fij) . Therefore, *E 2 (P) does not contribute to the damping rate 7 S (p); 
it will only shift the mass m s of the scalar. We thus can simply write: 

Is (P) = Im *£i Ko,p) + O (e 3 T) , (3.4) 

ZLU s q 

and henceforth, we will not manipulate *£ 2 (P) until we reach the calculation of the scalar 
energy uo s (p) in section four. 



A. Late momentum expansion 

In this subsection, we will first perform the Matsubara sums and analytically continue to 
real energies before expanding in powers of p. Using the structure ()2.4|) of the photon prop- 
agator in the strict Coulomb gauge and the expression of the one-photon-two-scalar vertex 
(|2.6|) . we see that *Ei(P) is composed of two terms. The first term involves longitudinal 
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photons and is denoted consequently by *Su(P), and the second one involves transverse 
photons and is denoted by *E lt (P). Let us look first at *£y(P). We have: 

*£ 1Z (P) = e 2 T £ y " [(2po - ^o) 2 (K) *A S (Q)] . (3.5) 



Writing explicitly the integral over the solid angle of k in a reference frame where p is the 
principal axis, we have: 

2 p+oo p+l 

*Z U (P) = —TJ2 / dkk 2 dx[(2p -k ) 2 *A l (k ,k)*A s (q ,q)}, (3.6) 



with x = k.p. Note that q = \J k 2 — 2pkx + p 2 , which makes the integral over x not feasible 
for the moment. 

We want to perform the Matsubara sum over k . For this, we use the spectral decompo- 
sition of the two dressed propagators. In general, we have [371 ]: 

pl/T p+oo 

Ai(k ,k)= dre lk ° T du Pi (u,k) (1 +n(w)) e - """, (3.7) 



where i stands for I, t or s and n(u) is the Bose-Einstein distribution. The explicit expressions 
of the spectral densities p { are displayed in (j3.12j) - (j3.15|) below. Making this replacement, 
the sum over /c can now be performed. The subsequent steps are standard: One imaginary- 
time integration is eliminated by a delta function and the second one yields an energy 
denominator. Everywhere except in the energy denominator po is replaced by 2nnT. The 
analytic continuation to real energies is taken at this stage and obtained by the replacement 
ipo — > uj s q(p) +i0 + . The imaginary part is extracted using the known relation 1/ (x + i0 + ) = 
Pr (1/x) — m5{x). We thus have: 

Im *£iz(P) = -^r~ [ dkk 2 [ dxf du Us0 — ^— p l (u, k) p s (uj s0 - u, q) . (3.8) 

In the above relation, only soft values of u are to contribute in the integral, and so we have 
used the approximation n (u) ~ T/ui. 

The transverse contribution *Si t is handled in similar steps and we obtain for its imagi- 
nary part the following expression: 

e 2 T r+°° r +1 ^ Q 
Im *T, lt (P) — p 2 dkk 2 dx (l— x 2 ) du — - — p t (uj,k) p s (u s0 —u,q) . (3.9) 

TT J v J -l 7-oo U(U S0 -LJ) 

Note that the transverse contribution to 7 S (p) already starts at order p 2 . 
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Now dividing by 2uo s q{j>) as required in (|3.4|) . we get the following longitudinal and 
transverse contributions to the damping rate: 



Isl (?) 
1st (P) 



' 8vr 



-OO y+1 

,2 



do; 



(2w s0 - loY 



pi (w, fc ) Ps Ko - w, g) 



2tt 



r; J — 1 ./ — oo 

dkk 2 l dx / duj \ \ 
1 J-oo W (W s0 - WJ 



Pt (w, fc) p s (wso - w, g) . (3. 10) 



According to (j3.4j) . the damping rate itself will be: 



7, (P) = Isl (P) + 1st (P) 



(3.11) 



Next we move to perform the integrals involved in the expressions of ()3.10j) above. The 
spectral densities p lt are known [12I 

p l>t (uj, k) = il>t (k) [5 (to - u lit (k)) -5(cu + u ht {k))\ + /3 M {u, k) 6 (k 2 - uj 2 ) , (3.12) 

where 3^ (k) are the residue functions and (3 l t (to, k) the cut functions. The residue functions 
are given by: 



u (uj 2 - A; 2 ^ 



li{k) 



k 2 (4/3 - uj 2 + k 2 



co (J 2 - k 2 ^ 



h(k) 



ui=uii(k) 



Auo 2 /3 - (uo 2 - k 2 Y 



(3.13) 



LL)=U>t{k) 



and the cut functions have the following expressions: 



/>,(*,«) = — 



2uo , fc + cj\ 2 47r 2 co> 2 



-1 



uj (k 2 - oo 2 ) 
3k 3 



k 2 - u 2 + 



2U 2 
3k 2 



k — to , k + to 
1 H In 

2kuo k — to 



+ 



2 , ,2\\2 



TYCO (k 2 — UJ 2 ) 



(3.14) 



Remember that all quantities are in units of m s . The spectral density p s does not have a 
cut. It simply writes: 



p s (uj, k) = 3 S (A:) [5 (to - u s0 (k)) - 5 (uo + ^ s0 (k))\ 



(3.15) 



with 3 S (A;) = l/2u; s0 (A;). 

We therefore replace the spectral densities by their respective expressions ()3.12|) and 
(|3.15|) . The integration over to disappears with the delta functions of (|3.15|) . From the ex- 
pression of p t t in (|3.12j) . there are going to be two kinds of contributions, i.e., a 5-contribution 
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and a G-contribution. Because of kinematics, the ^-contribution is always zero whereas the 
0-contribution survives. The integrations over x and k have to be performed numerically. 
One would fit the behaviors of (p) and 7 si (p) to several ultrasoft values of p, but in the 
spirit of the present work and in order to make direct comparison with the early momentum 
expansion method to be presented shortly, we display the results for j sl (p) and j st (p) in 
powers of p to fourth order. We obtain: 

e 2 T 

7 (p) = [1.44253 + 0.309278 p 2 - 0.133949p 4 + O (f)] + O (e 3 T) ; 

l07T L 



1st (?) 



e 2 T 

(1.65937 + 1.62114 In r?) p 2 

47T 



/ i Q4^37\ 

+ ( 71.3264 + 57.4152 In 77 + ' _ 2 j p 4 + (p 6 ) 



+ C(e 3 T). (3.16) 



£/m s and rj = i]/m s . 



The thermal mass m s has been reintroduced and here, p 

To order p 2 , these are the results obtained already in |36j]. Note that the longitudinal 
contribution 7^ (p) is safe from any infrared sensitivity whereas the transverse contribution 
7 st (p) is infrared sensitive, this despite the fact that both contributions have been handled 
in exactly similar ways. The reason why we have wanted to push the expansion to order p 4 
is to show that other forms of infrared sensitivity may appear, other than the familiar lnf/. 
Indeed, we clearly see the power-like behavior 1/fj 2 in the p 2 -coefficient of 7 st (p). These 
issues will be furthered in section five. According to (|3.11|) . the scalar damping rate is: 



e 2 T T 

7 On) = 1.44253 - (6.32820 + 6.48456 In r/)p 2 

167T |_ 

( 7 78148 \ 
285.172 + 229.661 In 77 + — — — J p 4 + O (p 6 ) 



+ C(e 3 T). (3.17) 



B. Early momentum expansion 



Now let us perform the same calculation while introducing from the start an early mo- 
mentum expansion, before the Matsubara sum and analytic continuation to real energies 
are done. First we perform analytically the integrals over the solid angle of k in the same 
reference frame where p is the principal axis. The difficulty comes from the presence of 
functions of q = \Jk 2 — 2pkx + p 2 , mainly the scalar dressed propagator *A S (qo,q). These 
need to be expanded and, in order to do this, we use the following expansion to fourth order 
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in the external momentum p: 



*A S (g , q) 



1 - px d k + 



V 



2 r 



1-x 2 



d k + x 2 d\ 



p 3 x 



1 - x' 
k 2 



1 ,fi 2 ™ 3 

-—r—dl - 

k k 3 fe 



4 

p 

+ 4fc 



-(l-6x 2 +5x 4 



2A: 



2fc 2 
+ (P 5 ) 



*A s (g ,fc). (3.18) 



Here stands for d/dk and, remember, x = k.p. Consider first *£ 1 ;(P). Insert the above 
expansion in its expression (|3.5|) and then perform the integrations over x which become 
straightforward. All odd orders in p cancel out and we are left with: 

/-+oo 



f^T J2 I dkk 2 (2p - hf * A, (k , k) 



1 + 



^ + ^)4(af + |a|) +0 ( P «) 



*A s {q ,k). (3.19) 



Only at this stage the Matsubara sum and analytic continuation to real energies are done. 
Here too we replace the dressed propagators by their spectral decompositions given in ()3.7|) 
and follow the standard steps sketched right after (|3.7j) in order to extract the imaginary 
part. Dividing by 2u s q (p) as indicated in ()2.12|) . we obtain the following expression: 



Isiip) 



' dkk 2 



4:71 



(2u s0 - u)" 



oo U LO (u s0 - to) 



Pi(u,k) 



x 



o{p 6 ) 



p s (uj s0 -uj,k). (3.20) 



The transverse-photon contribution *S lt (P) is handled in similar steps. Dividing by 
2uj s q (p), we obtain from it the following transverse contribution to the damping rate: 

du 

dkk" I — — 
/ /.){/. 

'v 



lst{P) 



O e 2rp f + OO f + OO 

I dkk 2 I 



LO (U s0 - u) 



p t {uj,k) 



x 



1 + 



2p 2 



p s (u s0 -u,k). 



(3.21) 



As one sees, there are different types of terms involved in ()3.20|) and ()3.21|) . In the sequel, 
we will briefly show how we carry through with each one. We will try to have the notation 
as clear and concise as possible. 
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1 . Integration 

The first type of integrals and the simplest we have to deal with is the one that involves 
a pp contribution with no derivatives. Generically, we consider an integral of the type: 



where i stands for I (longitudinal) or t (transverse). According to the form of 
the spectral functions (j3.12|) and ()3.15j) . there are two kinds of contributions: a 55 
(residue-residue) contribution and a Q5 (residue-cut) contribution. The 55 contribu- 
tion requires that the energies satisfy ±uj s0 (k) ± u>i(k) = 1, four constraints which 
are always forbidden by the dispersion relations. Hence the 55 contribution is al- 
ways zero because of kinematics. The Q5 contribution writes J^°°dk f^duf(k,ui,l — 
ui) 3s(/c)/3j (k, uj) [5 (1 — uj — uj s o) — 5 (1 — uj + oj s q)} 6 {k — \u\). A non-zero contribution 
must satisfy uj — 1 ± uj s q, together with —k < uj < k. It is not difficult to see that 
only the case uj — 1 — uj s o{k) is allowed, and this for all values k > i]. The integration over 
u is straightforward and we obtain for this type of integral: 

r+oa 

I PP = dki s f(k,l-u) s o,u) s o)/3 i (k,l-u s o). (3.23) 
Ji) 

Note that only integrals of the type l pp contribute to the coefficient of lowest order in p in 
the scalar damping rate. 

The second type of integrals we have to deal with is the following: 



Here too the discussion has to be carried out contribution by contribution, using the structure 
of the spectral functions (j3.12|) and (j3.15|) . The first contribution to consider is the one that 



When the derivative over k is applied to $ s {k), we have zero contribution. Indeed, in order 
for it to be nonvanishing, the energies must here too satisfy ±u s o(k) ± u)i(k) = 1, forbid- 
den by kinematics as already mentioned. But when applying the fc-derivative to the delta 
function, we also have zero contribution. This is simply because there is no intersection 
between the supports of the two delta functions, even if they are involved through deriva- 
tives, first order as in here or higher. The other contribution is a cut-residue, namely 




(3.22) 




(3.24) 
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Ir, dk I-oo duf{k,uj, 1 - uj)(5 i (k,uj)Q(k - \uj\) d k [i s 5 (1 - uj ^ u s0 )]. We first apply the 
derivative over k to j s and we get the piece dk},' s f{k, 1 — uj s0 , cj s0 ) Pi (k, 1 — uj s0 ), where 
l' s stands for d$ s (k)/dk. We then apply it to 5 (1 — uj =f ^so) an d use the standard rules regu- 
lating the handling of the delta distribution and we always check the results by regularizing 
either the derivative dk or the delta function itself. Only uo = 1 — uo s q contributes and we ob- 
tain f^°°dk uo' s0 3 S d u [f(k, 1 — uj, uj) (3^ (k, 1 — a;)] | w=a , s0 - Putting the above two contributions 
together, we obtain the following result: 

r+oo 

I P d kP = dk [$'J(k, 1 - ^ s0 , u s0 ) Pi (k, 1 - uo s0 ) 

+uj' s0h d u [f(k, 1-uj,uj) Pi (k, 1 - cj)}1 = „J . (3.25) 

The third type of integrals we have to deal with is one that involves a second derivative 
in k: 

r + oo p+oo 

I pd\p= I dk I dujf(k,uj,l-uj)p i (k,uj)dlp s (l-uj,k). (3.26) 

J T] J — OO 

The steps to treat the different contributions parallel those followed for I p d kP - As explained 
before, the SS contribution is zero because of kinematics. Therefore, the only contri- 
bution to look at is j^°°dk f^™dujf(k,uj, 1 — uj) Pi (k,uj) Q(k — \uj\) d\ \i s b (1 — uo =F uj s q)\. 
We have to use the identity d\ ($ s 5) = 3" 5 + 2$' s dk5 + 3 S <9f 5, where 3" stands 
for the second derivative of } s {k). The term involving i" s 8 yields J^ +oc 'dk$"f(k,l — 
^so, ^so) Pi (k, 1 — uj s0 ). The two other terms ^d^S + 3 S d%5 give together the following 
result j+°°dk [{ 1s uj" s0 + 2jXo) + h uof d 2 J [f(k, 1-uj,uj) P t (k, 1 - cj)} \ w=0Js0 . The proce- 
dure is to replace dk by uo S Qd u and apply the usual rules regulating the handling of the delta 
distribution. Putting the above results together, we obtain: 

r+oo 

I pd\p = dk [$"f(k, 1 - UJ s0 , Wso) Pi (k, 1 - uj s0 ) 

Jri 

+ [(3,^0 + 2 3 y s0 ) fc, + i s uj%dl] [f(k, 1-uj,uj) Pi (k, 1 - uj)} I ]. (3.27) 



Along similar lines we find the integral involving the third derivative in k: 

f + OO p + oo 

hdlp = dk duf(k, u,l- u;)pi{uj, k)dlp s (l - uj, k) 

J rj J — oo 

3 S 3) + (33'Xo + 3 3 X'o + 3^2) d u 
(33>f + 33^> s0 ) dl + Isu'lM [f(k, 1-uj,uj) pi (k, 1-uj)} }, (3.28) 



r+oo 
/ dk 
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and the one involving the fourth derivative: 



dk 



+00 



dk 



duf(k, uj,1- u)pi(u, k)dlp s {\ - uj, k) 
3i 4) + (43?^ + 63X0 + 43>S + 1*>§) 9, 



+ ( 6 3 y; + + h ( ^lu' s0 + 3j; ) ) &i + 1 4 3 y; Q + 63^^ ) 3 



+ h^sodu, \f(k,l-u,u)Pi(k,l-uj) 



\w=uisO 



(3.29) 



The next types of integrals we must handle are those coming from deriving the spectral 
density p s (uj s0 (p) — uj, k) in ()3.2()j) and (j3.21j) with respect to p. It is convenient to handle 
these terms by reintroducing 5 (oj s q (p) — uj — uj') and performing the derivative on it. The 
first such generic term is: 



dk 



duj 



dui' 8^5 (1 — uj — uj') f(k, uj, u^p^u, k)p s (u)', k) 



Tj J — OO J — OO 

hoo r+oo 

dk du d u [f(k,u,w')pi(k,u)]\^ =1 _^p a (l-u,k). (3.30) 

r\ J — oo 

Here too the residue-residue contribution is zero because of the same kinemat- 
ics. What remains to calculate is the cut-residue contribution, which is equal to 
- f^°°dk j^duiisdu [f(k,UJ,uj')(3 i (k,uj)Q(k - \u\)]\ u , = i_ u S (1 -w^w s0 ) and straightfor- 
wardly shown to yield: 



dki s d w [/(fc,cu,u; s o)A(^ w ) 



\U)=1— UJ S 



(3.31) 



The other term involves a second order u-derivative; it is worked out straightforwardly: 



l Ppd^S 



oo p+oo 

dk / du 

r\ J — oo •/ — -yL 

oo 



dJ d 2 J (l-oj- uj') f(k, u, uj^p^uj, k)p s (uj', k) 

dki s dl [/(fc,a>,u; a o)A(M)]L =1 _ w . • (3-32) 

Last are integrals involving an tu-derivative over 5(1 — uj — uj') and ^-derivatives over p s . 
The treatment is similar and straightforward too. With one fc-derivative we obtain: 

-+oo 

Ipdkpd^s — / dk 



Tj J — OO 

+oo 



duo I dw'd u 8(^ — u — uj')f{k,u,Lij')p i {uj,k)d^p s {uj',k) 



dk[t' s d u [f(k, uj, uj sQ )f3 i (k, u) 



Ioj=1— OJ* 



-lsu' sQ [d u [f(k,uj,uj s0 )p i (k,uj)} w=1 _ ujl ] | w/=w J 



(3.33) 
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and with a second-order fc-derivative we get: 

p+oo p+oo r+oo 

Ipdlpd^s = dk du dJ dj (1 - u - u/) f(k, u, uj')pi{uj, k)d%p s (u/, k) 

J Tj J — OO J — OO 

p + OO 

= - / dk[i" s ^ [f(k, u, UsQ^P^k, u)] | w= i_ WsQ 
-23X0 [/(^,^^ s o)/3i(^,^)] w= i_^] L, =Ws0 

-h Ko^ + [S w [/(*, w, ^ )A(fc, w)Uw] U J- (3-34) 

These are generically all the types of integrals we will have to evaluate. In each case, we 
have to examine the infrared behavior, i.e., the behavior close to rj. In the eventuality of the 
presence of an infrared divergence, we have to extract it analytically. The remaining finite 
part of the integral is generally performed numerically. 



2. Evaluation 



Let us work out a specific, general enough, example in some detail. Let us take it from 
the coefficient of p 4 in the expression in ()3.21|) for j st (p), namely the p t (u, k) dkp s (1 — uj, k) 
contribution. The corresponding function / is f{k,u,u)') = kj {uouo'). The integral involved 
is of the type I p d kP , see (I3.25J1 . Let us look at the first term $' s f(k, 1 — u s0 , u s0 ) (3 t (k, 1 — uj s q) 
where u> s0 stands here for u s0 (k). Using the expressions of /, the residue function $ s (k) 
given after f)3.15j) and (3 t (k, 1 — u s q) given in (|3.14|) . we perform a small- expansion of the 
product to obtain the following small-fc behavior: 

$'J(k, 1 - cu s0 , oj s0 )(3 t (k, 1 - lu s0 ) = -0.607927 jk + 4.70205 k + 30.0885 k 3 + O (k 5 ) . (3.35) 

It is clear that this is going to lead to a logarithmic divergence. To extract this divergence 
from the integral, we split the latter into two parts: one integration from 77 to any finite 
number plus one second integration from that finite number to +00. The integral writes 
then: 

p+oo 

/ dk$'J(k, 1 - u s0 , u s0 )(3 t (k, 1 - u s0 ) = 0.607927 In - 0.607927 lnk s0 (£) 

rh s0 (£) 

+ / dk [z'J{k, 1 - toso, to s0 )P t (k, 1 - u s0 ) + 0.607927/A;] 
Jo 

+ / dxi'J(k,l-u s0 ,uj s0 )P t (k,l-uj s o)\ k=kso{x y (3.36) 
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What we have done is this. In the finite part (the third term in (|3.36j) ). we have changed 
the integration variable from k to x = k/u s o(k), which implies that k = k s o(x). Note 
that k — > oo implies x — > 1. The finite value in question that splits the original inte- 
gral into two parts is then chosen in terms of x instead of k and is denoted by i. The 
second term in (|3.36|) is the original divergent piece of the total integral from the in- 
tegrand of which we have subtracted — 0.607927 jk in order to render it safe in the in- 
frared, and hence the lower bound rj is replaced by 0. We must then of course add the 
contribution -0.607927 f kso{£) f = 0.607927 In r] - 0.607927 In k s0 (£). As already men- 
tioned, the finite value t is arbitrary, but practically it must be chosen small enough in 
order to make the integral / ^ dk [$' s f(k, 1 — u s0 , u s0 ) (3 t (k, 1 — u s0 ) + 0.607927/A;] numer- 
ically feasible. Indeed, though we are assured of its finiteness analytically, both integrands 
$ r s f(k, 1 — u s q, Uso) (3 t (k, 1 — u s o) and 0.607927/fc are still each (here logarithmically) diver- 
gent for small k. However, when I is small, we can use a small-fc expansion in order to get a 
number for the integral. We have pushed the expansion to 0(k u ) and we start having good 
convergence for already i = 0.6. Also, we must (and do) check that the final result does not 
depend on a particular value of I. We finally get: 

/ dk$'J(k, 1 - u s0 , oj s0 )/3 t (k, 1 - oo s0 ) = 0.607927hi77 + 0.702549. (3.37) 

Jr) 

Now to the second contribution to the example we have chosen to detail, which is 
J^°° dk u}' s0 $ s duj [f(k, 1 — lj, lj) (3 i (k, 1 — ui)\ | w=Ws0 - A small-fc expansion of the integrand 
yields: 

<J*oh d * [f{k, I -to, to) p t (k, 1 - w)]L =a)s0 = -2.43171/fc 3 + 24.8687/A; - 208. 701 A; + 0(k 3 ). 

(3.38) 

Here we have a 1/r/ 2 behavior in addition to the familiar lnr/. We therefore write: 

r+oo 

/ dku' s0h d^ifik,! - (3 t (k,l - u)]l =UJs0 = -1.215 9/^ + 1. 215 

Jr, 

rk s0 (£) 

-24.86871n^ + 24.8687 In k s (£) + / dk [u' soh d„ [f(k, 1 - u,u)P t (k, 1 - u;)]^ 

Jo 

+ J dx u' s0 3 S [f(k, l-u,u) P t (k, 1 - Lj)]\ k=kso(x) . (3.39) 



2.43171 24.8687 



k 3 k 

We have proceeded as already explained. The integrations are smooth and no additional 
particular problem arises. Good convergence starts at £ = 0.6 and the ^-independence of 
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the sum is checked systematically. The final result is: 

r+oo 

/ dk cu' a0 h d u [f(k, l-u,u) /3 t (k, 1 - u)] \ U=UM = -1. 215 9/ V 2 - 24.8687 In 77 - 27.2104. 

Jri 

(3.40) 

Putting these two contributions together, we arrive for this specific example at: 

WLampic = -1-215 9/ry 2 - 24.26077 In 77 - 26.507851. (3.41) 

All the terms are treated along similar steps: infrared divergences are systematically 
detected by a small-/c expansion and extracted manually; the finite-part integrals performed 
numerically. We sum all the contributions and obtain exactly the two results ()3.16|) for 

Isi (P) and 1st (p), namely: 
e 2 T 

7 (p) = [1.44253 + 0.309278 f - 0.133949/5 4 + O if)} + O (e 3 T) ; 

167T L 



1st (?) 



e 2 T 

(1.65937+ 1.62114 In f]) f 

47T 



/ x 94537\ 

+ ( 71.3264 + 57.4152 In f] + — — — \ f + O (p 6 ) 



+ O (e 3 T) . (3.42) 



Note that by this method too the longitudinal contribution 7^ (p) is free from any infrared 
divergence whereas the transverse contribution is infrared sensitive. 



IV. SCALAR ENERGY 

Now we turn to the determination of the scalar energy u s (p) to next-to-leading order in 
the coupling constant e. It is defined in (|2.11jl . We want also to obtain it in both the late 
and early p-expansion. We will see that in this case too the same result is obtained, free 
from any infrared divergence. 



A. Late momentum expansion 

Recall that the next-to- leading order self-energy *E is the sum of two diagrams *Ej and 
*S 2 , see ()3.1|) . ()3.2|) and ()3.3|) . Recall also that we have stated that *S 2 is real. Indeed, using 
the structure of the two-photon-two-scalar vertex (|2.7j) and that of the photon propagator 
in the strict Coulomb gauge ([2.4)1 . we have: 

*£ 2 (P) = 2e 2 T^|^0[ *Ai(K) + 2 *A t (K)} . (4.1) 
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The angular integrals over the fc-solid angle and the sum over ko are done as explained 
before. No p-expansion is needed. We obtain the following expression: 

2rp r+co r+oo . 

*S 2 (P) = ± T dkk 2 — [ Pl (co, k) + 2p t (uj, k)} . (4.2) 

Jr, J-oo U 

As just restated, real. On the contrary, *Ex has an imaginary part which we have just 
manipulated in the previous section, both with late and early p-expansions. One useful 
way to obtain the real part of *£i is to use its imaginary part and the following dispersion 
relation K|: 

r+°° dt 1 

Re Tr * A M (k , k) *A S (q , q) f (k) = / [Im Tr * A M {ko, k) * A s (q , q) f (k)} . 

J-oo t-LOsQTl 

(4.3) 

Using this dispersion relation and performing the Matsubara sum before any momentum 
expansion, we obtain for the real part of the longitudinal contribution *Ey the following 
expression: 

Re*S 1 ,(P) = -^ [ + °dkk 2 f +C0 ^L_ ( +l dx [ + °° duJ ( 2u °° ~ u \ f Pl (uj,k)p s (t-0J,q). 



4vr 2 J v 7_ 00 t - uj s0 J_ x " u(t-u) 

(4.4) 

Remember, x = k.p. For the transverse contribution we obtain the following expres- 
sion: 

2rp P+OO p+OO i. n + 1 

Re *E lt (P) = -^v 2 \ dkk 2 — — / dx(l-x 2 ) 

7T 2 J v J-oo t ~ UJ s0 J -l 

/+oo ^ 
— -tp t (uj,k)p s (t-uj,q). (4.5) 
u{t-u) 

The next step is to perform the integrals in (|4.2jl . (|4.4jl and (|4.5jl . For *S 2 (P), we replace 
the photon dispersion relations by their expressions in ()3.12j) . The residue part leaves one 
integration over k and the cut part restricts the integration over u between —k and k. 
Otherwise the integrations themselves are done numerically. No early or late momentum 
expansion is necessary and no special problem arises. In particular, the infrared behavior is 
safe. Recall that we have to subtract the leading order terms for ultraviolet convergence as 
indicated after (|2.9|) . From the longitudinal photon we obtain — 0.183776 e 2 T and from the 
transverse photon — 0.000443967 e 2 T, a much smaller contribution. Summing the two, we 
get: 

*S 2 (P) = -0.184 22 e 2 T + 0(e 3 T) . (4.6) 
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No p-dependence as announced previously, which means *£2 (P) will only correct the scalar 
thermal mass m s , a correction free from any infrared divergence. 

To perform the integrations in (|4.4|) and (|4.5|) . we replace the scalar and photon spectral 
densities by their respective expressions ()3.15|) and ()3.12|) . Because the scalar spectral density 
does not involve a cut, the integration over t is trivial. Here too we could give values to 
the external momentum p and fit the resulting curves, but in the spirit of the present work 
and in view of the coming comparison with the results from the early momentum expansion 
method, we expand the integrand in powers of p and perform analytically the integration 
over x. The residue part in the photon spectral density will only leave the integration over 
k whereas the cut part will restrict the integration over to. The remaining single and double 
integrals are done numerically. Here too no special problem arises and the infrared region 
is safe. After subtracting the corresponding leading-order terms for ultraviolet convergence, 
we obtain: 

Re *E H (P) = e 2 T [-0.0746037 - 0.0201876 f + 0{f)\ + O (e 3 T) ; 

Re *E lt (P) = e 2 T [-0. 165908 f + 0(p 4 )] + O (e 3 T) . (4.7) 

It remains to sum (|4.6|) and (J4.7|) and divide by 2u) s0 (p) as instructed in (J2.ll)) to obtain 
the scalar energy: 

uj s (p) =m s [(1 + 0.258824 e + 0(e 2 )) + (0.5 + 0.056684 e + O (e 2 )) f + 0{f)\ . (4.8) 



B. Early momentum expansion 



Finally, let us recalculate the scalar energy u s (p) while introducing the expansion in p at 
an early stage, before performing the Matsubara sum. Since *S 2 does not depend on p, it 
will not be affected by this expansion; only *Y>i is affected. As in the calculation of 7 S (p), 
the integration over x becomes straightforward. We do the Matsubara sum and analytically 
continue to u s0 (p). We obtain for the longitudinal contribution: 

e 2 T /•+??. „ r +o ° dt r +co du 



Re*S i; (P) 



2tt 2 



dkk 2 



t — LU s0 



L)(t — U)) 



x t (2u sQ — uY 



p a (t-u,k). 



(4.9) 
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Odd powers in p cancel. The transverse contribution ()4.5j) is already of order p 2 and does not 
need further p-expansion. The x-integral done, Matsubara sum and analytical continuation 
to u) s0 (p) performed, we obtain: 

A 2rp r+oo p+oo i. p+oo j 

Re *Z lt (P) = ^p 2 / dkk 2 / / — f -t Pt (uj,k)p s (t-uj,k). (4.10) 

What remains now is to plug in the expressions of the spectral densities and perform the 
integrals over uo and t, and then over the momentum k. The scalar spectral density p s has 
no cut part, which makes the integration over t simple. The residue part in p l t eliminates 
the integration over u and the cut part restricts its limits. Numerical work is needed for the 
rest. Here follows a display of generic terms. First we have the type: 



oo r+oo p+oo 



R pp — I dk dt du>f(k,uj,t)p i (k,u!)p s (k,t — uj) 

J r\ J —oo J —oo 

i 

dk \$ifa[f(k,Ui,Ui + u s0 ) - f(k,Ui,Ui - u) s0 ) - f(k, -Ui, -Ui + u s0 ) 
~+k 



v 



/+« 
du^fa (k,uj) [f(k,u,u + uj s0 ) - f(k,uj,uj - cj s0 )]].(4.11) 
■k 

The subscript i stands for I or t. Note that only integrals of type I pp contribute to the 
coefficient of zeroth order in p in the real part. The second type of integrals is the following: 

r+oo r+oo f — oc 

R pdkP = dk dt dujf(k,u,t) p i {k,u)d k p s {t -uj,k) 

J t] J —oo J —oo 

" + 00 

dk [iii' s [f(k,Ui,Ui + uj s0 ) - f(k,Ui,Ui - uj s0 ) - f(k, -u h + uj s0 ) 



V 

+ f{k, -uJi, -uji - u s0 )] + u' s0 ^i s [d t UXk,Ui,t) + f{k, -Ui, -t))\ t= ^ +uJso 

-k 

-f(k,u,u- u s0 )} + J s0 tsPi (k, u) [d t f(k, w, t) \ Mujs0 + d t f(k, u, t) \ t= ^J}. (4. 12) 
The third and last type of integrals is the one involving a second-order derivative in k; it 
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writes: 

-+oo r+oo r+oo 



r+oo r+oo r+oo 

R P d 2 kP = dk dt dwf(k,u,t)p i (k,u)dlp s (t-uj,k) 

J n J — oo J — oo 

dk [ h fc' + u% is d 2 t ) [{f{k,u h t))\ t= ^ M - (/(fc.Wi.t))!^ 



'/ 



x [d t (f(k, Ui,t) + f(k, -u u -t))\ t= ^ s0 + d t (f(k, u,i,t) + f(k, -ui, t))\ t=UJi _ w J 

/+k 
-k 

+ (a^ + o&jjft (k,u) [9tf(k,uj,t)\ t=li)+Ua0 + d t f(k,u,t)\ t= ^J]] . (4.13) 

Using these generic results, we perform the integrals numerically. We do not encounter 
any additional problem and we find all the integrals safe in the infrared. Here too we have to 
subtract the leading-order terms for ultraviolet convergence and, putting all contributions 
together, the numerical integrations yield the result: 

uj s (p) =m s [(l + 0.258824 e + 0(e 2 )) + (0.5 + 0.056684 e + O (e 2 )) f + 0(p 4 )] . (4.14) 

This is exactly the same result (|4.8jl obtained with the late momentum expansion. We should 
note the fact that there is no infrared sensitivity in the above expression, even though it has 
been technically more involved to derive than the determination of the damping rate 7 S (p). 

V. DISCUSSION 

This work aimed at calculating the damping rate 7 S (p) and energy u s (p) for scalars in the 
context of next-to-leading order hard-thermal-loop perturbation theory for scalar QED. For 
each of the two quantities, we have carried out the calculation in two ways: using both a late 
and early external momentum expansion, respectively after and before the Matsubara sum is 
done and the analytic continuation to real energies taken. Though technically different, both 
methods yield the same results for the damping rate and energy. This fact is particularly 



interesting in view of an early work [33J, |34j in which the early momentum expansion was used 
in the context of next-to-leading order hard-thermal-loop perturbation of hot QCD in order 
to extract analytically a value for the non-moving longitudinal-gluon damping rate. The 
very definition of that quantity imposed a systematic expansion in powers of the external 
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momentum p, at least to second order. From a technical point of view, the question was 
where to perform the expansion. Ideally, the latest possible, at least after the Matsubara 
sum and the analytic continuation to real energies are done. But technically that was 
unfortunately not feasible in QCD because of the complication of the intermediary steps, at 
least if one wanted to carry through analytically. The early momentum expansion was used 
instead, and legitimately the occurrence of infrared sensitivity of the result was partially 
incriminated on the method used. 

The present work shows that infrared divergences can occur without the early momentum 
expansion. This is clear in the transverse contribution 7 st (p) to the scalar damping rate, see 
(|3.16|) . Remember that up to order p 2 , this result with the logarithmic divergence is already 
found in [36j] . The reason we have pushed the calculation to order p 4 is to demonstrate that 
other forms of divergences do occur, here a 1/rf. This was the case in QCD too. Since 
these results are obtained without the early momentum expansion and the same results are 
obtained with it, we think we have here a strong indication that this latter is not responsible 
for any infrared sensitivity. 

Furthermore, infrared sensitivity is not systematic when we perform a momentum expan- 
sion, late or early. This is exemplified in the longitudinal contribution 7 s/ (p) to the damping 
rate, see f)3. . and, maybe more pertinently since the calculations are more intricate, in the 
scalar energy uj s (p). It is true that only fourth order in p is included in the determination 
of 7 S ; (j>) and second order in the determination of uj s (p), but the trend of the calculations 
indicates that higher-order coefficients will eventually be safe. 

A couple of objections may still rise. First, one could argue that in fact, the late and early 
momentum expansions are actually the same for the quantities we have treated since the 
expressions of these, before any momentum expansion is performed, are still not analytically 
closed in terms of p . To answer this objection, one can go to the photonic sector of scalar 
QED treated in Indeed, the HTL-summed next-to-leading order photon self-energies 
are evaluated in closed form for all ui and p. Let us focus only on the longitudinal photons. 
Three regions in u and p are to be distinguished: u 2 < p 2 (region I), p 2 < uo 2 < Am 2 + p 2 
(region II) and Am 2 + p 2 < u 2 (region III). The longitudinal next-to- leading order HTL- 
summed photon self-energy is found to have the following expression 



e 2 T u 2 - p 
8n p 2 



,2 r 



4m s + 2ie — i — In 
p 



( 



2m s — i (e — p) 



2m s — i (e + p) 



) 



(5.1) 
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where the prefix 8 indicates here too that the known leading hard-thermal-loop contribution 
has been subtracted for ultraviolet convergence, and e = Q in regions I and III and e — i |G| 
in region II, with G 2 (u),p) = u 2 (u 2 — p 2 — Am 2 ) / (u 2 — p 2 ). The full longitudinal-photon 
dispersion relation is: 

n 2 = p 2 + u l (n hP ), (5.2) 

where 11/ is the full longitudinal-photon self-energy, and, up to next-to-leading order, writes: 



&l IP) = P + ~, o rri 7 V' 5 ' 3 

where 511/ (u,p) is the hard thermal loop given by [l3j: 

5U t ( w , p) = 3ml (l ~ ^) (l - 1 1" ) . (5-4) 

with m p = eT/3 the photon thermal mass and cjj (p) the on-shell longitudinal photon energy, 
solution to ()2.8|) to lowest order where only the hard thermal loop is kept in the self-energy. 
Using ([5.4)1 . we can rewrite (|2.9J) as: 

n? (p) = p 2 + ^/ + - ' 2 lPJ - rn ; ( UhP ) . (5.5) 

d?7lp + p z — UJJ (p) 

Remember that all these results involve no expansion in p and the Matsubara sum and 
analytic continuation to real energies are already done. Also, all intermediary integrals are 
performed. Now we expand. The region of interest to us is region II, where we are allowed 
to perform the following expansion: 

3 3 

uji (p) = 1 + — p 2 - — p A + O (p 6 ) ; p = p/m p . (5.6) 

Using this, we perform the expansion of ()2.10|) in powers of p. It is straightforward and we 
find: 



n 2 (p) = m 2 p 



(1 -0.368e + O(e 2 )) + - 0.0536 e + O (e 2 )^j p 2 + 0(p 4 ) 



(5.7) 



Here we have together the leading and next-to-leading orders in the coupling e. Note that 
there is no infrared sensitivity. What remains now is to perform the expansion in powers of 
p before the Matsubara sum and analytic continuation to real energies are done, and carry 
on with steps similar to what we have gone through in this work. This we do and we obtain 
exactly the same result (|5.7jl . There is no need to display the results we think. This is an 
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additional indication that the early momentum expansion works. Of course, when it can 
be avoided like in the present context of scalar QED, there is no need to use it. But in 
the context of hot QCD, a more realistic theory, it is sometimes necessary if one wants to 
manipulate analytically. 

Last, one may argue that if the calculations are to lead to infrared sensitivity, then the 
regularization with an infrared cut-off rj should not be used in the first place. To answer this, 
one can say that, first of all, infrared sensitivity is not systematic in all similar quantities that 
are treated in similar methods. We have seen a good example of this in this work, which is 
the two contributions to the scalar damping rate 7 S (p): the longitudinal contribution j sl (p) 
is infrared safe and the transverse contribution 7 si (p) infrared sensitive, though both are 
calculated in exactly the same manner. Hence infrared sensitivity may be a feature of the 
physical quantity in question. Furthermore, we could have taken the transverse contribution 
7 st (p) , made no expansion in p and tried to determine it numerically for different values of p 
without introducing any infrared cutoff. Then the internal momentum /c-integration would 
not converge, because precisely of the very infrared sensitivity. One could argue that the 
regularization should be different from the naive introduction of an infrared cut-off. But, in 
whatever way it is chosen, regularization will only exhibit a potential infrared sensitivity; it 
will not remove or cancel it. For that purpose, some other procedure has to be invoked, and 
we have never meant to do that. This particular issue would be pursued elsewhere. 
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